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Abstract. Hard partons propagating through hot and dense matter loose energy, leading to the observed 
depletion of hard hadron spectra in nucleus nucleus collision as compared to scaled proton proton collisions. 
This lost energy has to reappear in the medium due to the conservation of energy. Apart from heating the 
medium, there is the possibility that a propagating collective mode is excited. We outline a formalism that 
can be used to track the propagation of such a mode through the evolving medium if its dispersion relation 
is known and to calculate the resulting distortion of hadronic spectra at freeze-out. Using this formalism, 
we demonstrate within a detailed picture of the evolution of the expanding system and the experimental 
trigger conditions that the assumption of a sound mode being excited is in line with 2-particle correlation 
measurements by PHENIX and STAR for a semi-hard trigger. In this case, the measurement is sensitive 
to the averaged speed of sound in the hot matter. We also make suggestions how this sensitivity can be 
improved. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 25.75.Gz Particle correlations 



1 Introduction 

Energy loss of a high pt 'hard' parton travelling through 
low pt 'soft' matter has long been recognized as a promis- 
ing tool to study the initial high-density phases of ultrarel- 
ativistic hevay-ion collisions (URHIC) 1,2,3,4,5,6 . How- 
ever, if one considers the whole dynamical system created 
in the collision of two relativistic nuclei and not only the 
partons emerging from a particular hard scattering ver- 
tex, energy is not lost but rather redistributed into the 
medium. 

Measurements of angular correlations of hadrons as- 
sociated with a given hard trigger allow, as a function of 
associate hadron momentum, to study how and at what 
scales this redistribution of energy and momentum takes 
place. Such measurements for semi-hard hadrons with 1 
GeV < pt < 2.5 GeV associated with a trigger 2.5 GeV 
< Pt < 4.0 GeV have shown a surprising splitting of the 
away side peak for all centralities but peripheral collisions, 
qualitatively very different from a broadened away side 
peak observed in p-p or d-Au collisions For a harder 
trigger 2.0 GeV < pt < 4 GeV and different associate 
hadron momentum cuts, a persistent large angle signal 
has been observed [B]. 

As most promising explanation for these findings, the 
assumption that Mach Shockwaves are excited by the en- 
ergy lost from the hard parton to the medium has been 
brought forward [MTU] . In the following, we investigate un- 
der what conditions signals from such Shockwaves remain 
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observable in the dynamical environment of a heavy-ion 
collision, provided that a realistic description of the evolv- 
ing medium and the experimental trigger conditions are 
taken into account. Here, we summarize and expand on 
Q31 and [12]. 

2 Shockwave excitation, propagation and 
freeze-out 

For the time being we focus on central collisions only. 
For the description of the evolving medium, we employ 
a parametrized evolution model |13U14j which is known to 
describe bulk matter spectra and HBT correlations well. 
Since we're interested in the deposition of lost jet energy 
into the medium, our first task is to determine the space- 
time pattern of energy loss. We use the quenching weights 
from rjS] to obtain the probability P(AE) to lose the 
amount of energy AE from the two key quantities plasma 
frequency 

w c (ro,0 = f mm (!) 
Jo 

and averaged momentum transfer 

(qL)(r ,<f>) = [ d£q{£) (2) 
Jo 

in a static equivalent scenario which are calculated 
along the path of the hard parton through the medium. 



2 



Thorsten Renk: Signatures of Mach shocks at RHIC 



Since we are not interested in folding the result with a 
steeply falling spectrum but rather into the energy de- 
posited on average in a given volume element we focus 
on the average energy loss (AE) = J °° P(AE)AEdAE in 
the following. 

We make no attempt to calculate the microscopic mech- 
anism by which a propagating mode is excited. In [Hj it was 
shown that only a certain class of source terms leads to a 
propagating sound wave in a hydrodynamical medium. In 
particular source terms which directly couple energy and 
momentum lost by the hard parton into the local fluid cell 
do not lead to angular correlations as seen by experiment. 
This is in fact consistent with the findings of where no 
such correlations were found in a 2d Bjorken hydro code 
under the assumption of a direct coupling. A possible mi- 
croscopical excitation mechanism might be colour wakes 
induced in the medium |17| . 

For the purpose of the present work, we assume that 
a fraction / of the energy and momentum lost to the 
medium excites a Shockwave characterized by a dispersion 
relation E = c 2 p and a fraction (1 — /) in essence heats 
the medium and leads to some amount of collective drift 
along the jet axis to conserve longitudinal momentum. 

We calculate the speed of sound c s locally from a quasi- 
particle description of the equation of state as measured 
on the lattice [HJ as c s = dp(T)/de(T). This EOS shows 
a significant reduction of c s as one approaches the phase 
transition but doesn't lead to a mixed phase. The dis- 
persion relation along with the energy and momentum 
deposition determines the initial angle of propagation of 
the shock front with the jet axis (the 'Mach angle') as 
4> = arccosc,,. We discretize the time into small inter- 
vals At, calculate the energy deposited in that time as 
E(t) = At ■ dE/dr we then propagate the part of the 
shockfront remaining in the midrapidity slice (i.e. in the 
detector acceptance) . Each piece of the front is propagated 
with the local speed of sound and the angle of propagation 
is constantly corrected as 



matching conditions, we use the Cooper-Frye formula 



arccos ■ 



fc E C a {r)dT 
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where c s (r) is determined by the propagation path. 

Once an element of the wavefront reaches the freeze- 
out condition T — Tp, a hydrodynamical mode cannot 
propagate further. We assume that this point that the 
energy contained in the Shockwave is not used to produce 
hadrons but rather is converted into kinetic energy. In the 
local restframe, we then have a matching condition for the 
dispersion relations 



E = c 2 s p and E = y/M 2 



M 



(4) 



where M — V (p(Tp) + e(Tp)) is the 'mass' of a volume 
element at freeze-out temperature. 
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to convert the fluid element into a hadronic distribution. 
The resulting momentum spectrum is thus a thermal two 
component spectrum resulting from an integrations in- 
volving volume not part of the Shockwave and volume 
receiving an additional boost from the Shockwave. 



3 Simulation of the trigger conditions 

We simulate the PHENIX trigger conditions as closely as 
possible using a Monte-Carlo approach. We start by gener- 
ating vertices with a distribution weighted by the nuclear 
overlap 



T AA (h) = J dzp 2 (b,z) 



(6) 



a volume element receives from the Shockwave using the 



We then determine the jet momentum and parton type 
by randomly sampling partonic transverse momentum spec- 
tra generated by the VNI/BMS parton cascade as de- 
scribed in rjHj- Calculating the energy loss of the near 
side parton, we decide if the experimental trigger condi- 
tion is fulfilled. Since the experiment triggers on a hard 
hadron in the transition between the recombination and 
fragmentation regime, the model at this point cannot im- 
plement the trigger condition exactly. Instead, we require 
the trigger condition to be fulfilled by the parton and have 
checked that the model results do not change significantly 
when the near side trigger threshold is increased by 2 GeV. 
We note that this procedure places the vertices fulfilling 
the trigger condition close to the surface of the produced 
matter, i.e. in our model the medium is rather opaque, in 
agreement with the conclusions of [20ll21| . 

Once the vertex and momentum of a near side jet has 
passed the trigger condition, we determine the direction 
of the away side parton in the transverse plane and rapid- 
ity. In order to take into account intrinsic kx, we do not 
propagate the away side directly opposite to the near side 
jet but allow for a random angle. We have verified that 
this distribution, folded with the width of the near side 
peak reproduces the width of the far side peak in the case 
of d-Au and 60-90% peripheral Au-Au collisions. 

With vertex, energy and direction of the away side jet 
fixed, we calculate dE/dr of the outgoing parton. We stop 
the calculation when a significant fraction of the energy is 
lost to the medium. In each event we assume that a frac- 
tion (1 — /) of the energy lost from the hard parton heats 
the medium. We assume that due to momentum conserva- 
tion this contribution not leading to Shockwaves will lead 
to additional flow into the direction of the original away 
side parton. Likewise, we account for the possibility of a 
punchthrough if the initial vertex is very peripheral and 
both near and away side parton propagate near-tangential 
to the surface. 
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4 The role of transverse flow 

It has been argued that the includion of transverse flow 
leads to a substantial distortion of the Mach cones [22] ■ We 
observe this effect when we consider the spatial position of 
the cone at a given time However, as apparent from 
Eq. (|5J), the measured effect of the Shockwave is not man- 
ifest in position space but rather (via a boost to volume 
elements undergoing freeze-out) in momentum space. 

When the associate hadron cut is set well above the 
bulk matter momentum scales at freeze-out p rs 3Tp ~ 
400 MeV, the observation takes place in the tail of the 
boosted thermal distribution. Under these conditions, the 
yield is very sensitive to boosts and almost all observed 
yield comes from a region where the transverse flow is 
aligned with the emitted particle momentum. In the pres- 
ence of a Shockwave, the signal at 1 GeV momentum is 
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jiow are p ara ji e i (j n f a ct, the ob- 
'"" '' 1 "if} ow is about 9 times larger 
than the yield for configurations in which u s ^ ock _L ) . 
However, if flow velocity and shock propagation are (ap- 
proximately) aligned, no distortion of the Mach angle can 
take place. Thus, for an associate momentum cut beyond 
typical thermal scales, only configurations in which the 
flow does not distort the Mach angle are visible. 

On the other hand, momentum conservation dictates 
that the correlation signal cannot simply vanish for con- 
figurations in which shock and flow are not aligned. In this 
case, a broader peak structure is recovered with a lower 
Pt cut. This is illustrated in Fig. Q 



5 The role of longitudinal flow 

Even if the trigger parton is confined to be at midrapidity, 
the rapidity of the away side parton is not known. The 
differential production cross-section of hard partons in A- 
A collisions can be obtained in leading order pQCD by 
folding the two particle production cross-sections da/dt 
with the nuclear parton distribution functions fij/A (here 
we use 1221: Q 2 dependence has been suppressed in the 
following expressions): 



( p (J AA-*kl+X 

dp^dyidy 2 



d&v 



h3 



dt 



(7) 



If the outgoing partons are at rapidities y\ and y 2 , x\ 
and x 2 are determined by: 



Xl,2 



-j= [exp(±yi) 
\/s 



exp(±y 2 )] 



The conditional probability distributions P(y) of pro- 
ducing an away-side parton at rapidity y can then be cal- 
culated from the normalized cross-section (Eq. (7J) given 
the trigger parton at y\ — 0. For the dominant processes, 
it is a rather wide plateau between — 2 < y < 2. 

This raises an obvious question: Since the Mach cone is 
symmetric around the away side parton direction, a mea- 
sured large angle correlation near midrapidity also implies 



large rapidity correlations at zero transverse angle. Thus, 
the averaging over P(y) will tend to smear the signal mea- 
sured at midrapidity out towards smaller angles as com- 
pared to a simple midrapidity projection of teh cone for 
an away side parton emerging at y = 0. 

However, since the shock wave travels with c s in the 
local rest frame, the spatial position the of the shock front 
has to be determined by solving the characteristic equa- 
tion: 



dz 
It 



u{z,R, t) + c s (T(z,R,T)) 



z=z(t) 



l + u(z,R,t)c s (T(z,R,t)) 



(9) 



z=z(t) 



We use this equation to infer the longitudinal boost 
for an element of the Mach cone at forward rapidity. In 
essence, this means that a Mach cone in tfi, y-space is 
elongated significantly in y direction by longitudinal flow 
because a sound wave (unlike emitted particles) propa- 
gates in rapidity space even in a Bjorken expansion. For a 
Bjorkcn expansion this elongation amounts to a about 1.5 
units rapidity. The net effect of this elongation is that the 
Mach signal, i.e. the angular correlation at given rapidity 
are much less sensitive to the rapidity averaging procedure 
than naively expected. We illustrate this in Fig. 
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Fig. 2. Calculated 2-particle correlation under the assumption 
that a) the away side parton is always at midrapidity and the 
excited mode doesn't couple to flow (green) b) using realistic 
P{y) and assuming that the excited mode doesn't couple to 
flow (red) and c) including realistic P(y) and longitudinal flow 
elongation. 



(8) 6 The dip 



It is sometimes argued that it is the dip at the expected po- 
sition of the away side parton (at degrees in Fig.[5J along 
with the rise of the correlation strength at some large an- 
gle with respect to this which constitutes the signal of 
Mach shocks. However, even in our analysis, the fraction 
of energy exciting the Shockwave is only 75%, thus 25% of 
the energy lost by the away side parton ends up associated 
with its original direction. It is easily apparent from Fig. 3 
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Fig. 1. The effect of transverse flow on the observable correlation peak for configurations in which transverse flow and cone are 
aligned (solid) and orthogonal (dashed, 6 = 3 fm) for associate hadron pr = 1.0 GeV (left panel) and 500 MeV (right panel). 
The apparent disappearance of the correlation signal at high pr when flow and Shockwave are not aligned is clearly visible. 
At lower pt the signal reappears significantly broadened, reflecting the distortion of the Mach angle by transverse flow. All 
calculations here are for a single event with given parton energy E and assuming no intrinsic kr and y = as the position of 
the away side parton. Correlation strength in the direction of the away side parton have been suppressed for clarity. 



in JTj that the appearance of the dip is strongly linked to 
the fact that the dominant fraction of lost energy excites a 
Shockwave, not to the fact that a Shockwave is present as 
such. It is equally instructive to study Fig.Elin the present 
work. The same signal that would be interpreted as having 
a dip if the rapidity averaging is carried out fully would 
change dramatically if the position of the away side parton 
could be confined to midrapidity. The reason is that while 
in the rapidity-averaged situation the away side parton is 
outside the narrow acceptance region around midrapidity 
in most events (and hence only a part of the wider ring 
ends up within the acceptance), for specified away side 
parton rapidity suddenly all correlation strength associ- 
ated with the original direction ends up always inside the 
acceptance region whereas the Mach ring is still for the 
most part outside. 

Thus, the presence or absence not only hinges on the 
question if the cone is the dominant mode excited by lost 
energy but also on the experimental cuts. It is also ex- 
pected that for increasing trigger energy a larger fraction 
of original away side partons emerges from the medium, 
subsequently undergoing vacuum fragmentation. Such a 
contribution would likewise fill the dip region. 

A more reliable characteristic of the Mach cone is that 
the angle to which the away side correlation extends is 
a property of the medium (it is solely determined by the 
path-averaged speed of sound, cf. Eq. |J3J and independent 
of both trigger energy and associate hadron cut. Such a 
behaviour is difficult to obtain in radiative processes from 
the away side parton where it is expected that for increas- 
ing kx of emitted secondaries (and hence increasing asso- 
ciate hadron cut) the angle shrinks. A signal qualitative 
in agreement with this expectation of unchanged angle is 
observed by the STAR collaboration [Hj. 



7 Precision determination of the Mach angle 

If the observed correlation pattern is indeed caused by 
Mach shocks, a determination of the precise Mach angle 
offers an opportunity to probe the equation of state of 
the underlying system. In |1 1 j we have shown that there 
is some sensitivity to the averaged speed of sound - the 
difference between assuming a phase transition with a soft 
point (with small c s ) and assuming c s w 0.55 at all times 
amounts to a shift of some 18 degrees in peak position. 
However, more realistic variations of (c s ) accessible by ex- 
periment (e.g. by a variation of initial temperature for 
different y/s or a variation of decoupling temperature by 
going to different system size) to confirm our ideas about 
the medium EOS require a determination of the peak posi- 
tion by a precision of the order of a few degrees (0(1 — 5°)). 
This is beyond the ability of the current data and made 
difficult by the fact that in the present measurement the 
peak is always broadened. 

In a static, homogeneous cold medium, the correlation 
signal would closely resemble a (S-function at 4>Mach (which 
is also different from 4>Mach in the dynamical evolution). 
In the dynamical model, the correlation signal appears 
broadened for a number of reasons: 1) Different elements 
of the cone are excited at different times and probe differ- 
ent parts of the (inhomogeneous) medium. This is disper- 
sion introduced by the spacctime structure of the medium. 
2) The intrinsic kx randomly influences the angular po- 
sition of the away side parton. Thus, in averaging over 
many events all peaks get smeared with the induced an- 
gular spread. 3) The thermal motion of hadrons at freeze- 
out leads to a thermal smearing of the correlation signal 
and 4) The rapidity averaging outlined in section [5] leads 
to a widening and systematic shift towards small angles 
when averaging over many events is done. 

However, some of these effects can be addressed by 
experiment: By requiring a hard dihadron as the trigger, 



Thorsten Renk: Signatures of Mach shocks at RHIC 



5 



i.e. a near side hadron and a (softer) away side hadron, 
the rapidity and intrinsic fc^-kick of the away side par- 
ton is known for each event and the Mach correlation can 
be studied relative to the measured away side parton di- 
rection. In doing so, 2) and 4) are eliminated as mech- 
anisms leading to a widening of the peak. The thermal 
width can be somewhat reduced by increasing the mo- 
mentum threshold of analyzed associate hadrons (thermal 
motion transverse to the associate hadron direction is of 
less relative importance if the associate hadron momentum 
is increased). However, one can not push this to arbitrar- 
ily high momenta as the presence of a Shockwave signal 
requires that a fluid medium dominates the hadronic spec- 
tra (which ceases to be true above 2-3 GeV). The resulting 
peak in the correlation signal could be significantly nar- 
rower than what is observed now (cf. Fig. ^ where intrinsic 
kx broadening and rapidity averaging are absent since the 
plot shows a single event only). 

Finally, measuring correlations for hadrons with a dif- 
ferent freeze-out systematics (e.g. <j> an d £2 for which blast 
wave fits and their small scattering cross suggest early de- 
coupling) could potentially offer a window to study both 
the systematic change in thermal width of the signal and 
a different Mach angle (since Eq. ® is evaluated for a 
different upper limit). 



8 Conclusions 

While we have not presented a complete microscopical 
theory of the excitation of Mach Shockwaves by hard par- 
tons travelling through a hot and dense medium, we have 
shown that if one makes the assumption that such shock- 
waves are excited, the resulting correlation pattern is con- 
sistent with the measured data even if a realistic trigger 
simulation and effects of the medium evolution are taken 
into account. 

In particular, we have presented the (somewhat coun- 
terintuitive) finding that above a certain threshold in trans- 
verse momentum, the correlation signal is not significantly 
distorted by transverse flow. The underlying reason is pe- 
culiar for a hydro phenomenon — since scales 0(1 GeV) 
where the correlation has been measured by PHENIX are 
above the typical hydro scales, only alignment of flow and 
Shockwave direction (and hence no distortion) can boost 
the hydro medium enough to create a signal. 

We presented a similar finding within a discussion of 
the longitudinal dynamics. There, the elongation of the 
Mach cone by longitudinal flow, a phenomenon only tak- 
ing place for a perturbation moving with a given speed 
relative to the expanding medium, proved to be a crucial 
ingredient in the survival of the signal after averaging over 
the (a priori unknown) position of the away side parton 
in rapidity. 

While a detailed discussion of the excitation function 
as a function of trigger energy and associate cut still has 
not been done (the chief obstacle being a clean modelling 
of the transition between recombination and fragmenta- 
tion as a mechanism for hadronization) we note that the 



absence of cone-like structures in high transverse momen- 
tum angular correlation measurements |24U25| and their 
explanation in terms of hard partons emerging from the 
medium |21| are well in line with the observation that 
hydrodynamical modes are irrelevant at px > 4 GeV at 
RHIC. 

In summary, while the current data certainly do not 
allow Mach shocks as the only explanation for the mea- 
sured correlation pattern, it can be shown that they are 
a consistent explanation, and their nature as a hydrody- 
namical mode explains some characteristics found in the 
data in a natural way. A similarly conclusive case has yet 
to be made for alternative explantions for the correlation 
pattern seen in the data. 
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